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Q \ Abstract 

^ ' The 2/5 decay study of ^^'^Gd has been performed in the Solotvina Underground Laboratory 
^ . with the help of Gd2Si05:Ce crystal scintillator (volume 95 cm^). The background of the 
detector in the vicinity of the energy of ^^°Gd was reduced to 1.0 cpd/keV-kg. The new 
^ \ improved half-life limits have been established for Oz/2/5 decay of ^^°Gd to the ground (0^) and 
O ; first excited (2+) levels of ^^^Dy: TO/2 > 2.3(1.3) xlO^i yr at 68%(90%) C.L. The T1/2 bounds 
^ ' have been also set for 2i/, Oz/x and Oz/^X modes of ^^°Gd decay, as well as for different 2/? decay 
processes in ^^^Ce, ^^'^Ce and ^"'^Ce. 
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O : PACS: 23.40.-s; 12.60.-i; 27.60.+j 
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C 1 Introduction 

\ An exceptional interest to the neutrinoless double beta (Oz/2/5) decay is explained by the great 
\ potential of this process - which violates the lepton number conservation - to search for the 
neutrino mass {mi) and its nature as a sign of a possible new physics beyond the standard 
model (SM) [Q, ^ H, ^, |^, |^. The absence of the Oi/2/3 decay, established at the required level 
of sensitivity, yields strong restrictions on m^,, lepton violation constants and other parameters 
of the manifold SM extensions, which allow to narrow a wide choice of theoretical models and 
to touch the multi-TeV energy range competitive to the accelerator experiments [^, ^, ^, ^. 

In fact, the most sensitive Qv2(3 results were obtained by using the so called "active source" 
technique, in which a detector (containing 2/5 candidate nuclei) serves as source and detector 
simultaneously [|l], 0]. As examples, we recall the impressive half-life limits T^j2 in the range of 
(1-4) X 10^3 yr (m^ < 1.3-4 eV) established for ^^^Xe (high pressure Xe TPC) §, ^^^Te (low 
temperature bolometers Te02) pl, ^^^Cd (enriched ^^^CdW04 scintillators) 0, and the highest 
limit reached for ^'^Ge (enriched HP ^'^Ge detectors): T^J^ > 1.8 x 10^^ yr 0, Therefore, 



^Corresponding author. Institute for Nuclear Research, Prospekt Nauki 47, MSP 03680 Kiev, Ukraine; fax: 
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it is apparent that application of the "active source" technique with a new nucleus allows to 
extend the number of 2(3 candidates studied with a high sensitivity^. 

During last years cerium-doped gadolinium silicate Gd2Si05:Ce (GSO) crystal scintillators 
have been developed |T3|, |T^ . These scintillators are non- hygroscopic and have a large density 
(6.71 g/cm^), fast response (primary decay time about 30 — 60 ns), quite high light output (20% 
of Nal(Tl), wavelength of emission maximum 440 nm). Moreover, it was already demonstrated 
|T5|, |16], [1^ that GSO crystals can be applied for the 2/? decay search of ^^°Gd, which is one 
of the interesting candidate nucleus due to the following reasons. First, despite rather low 2/3 
decay energy release {Qpf^ = 1729.7(13) keV pl) its theoretical value ofT^)'2-{m^f = 8.6 x 10^^ 
yr is nearly three times lower than that for ^^Ge and ^^^Xe |T^, thus for the equal measured 
limits the experiment with ^^°Gd will yield the more stringent restrictions on the neutrino 
mass and other parameters of the theory. Secondly, recent calculation shows that two- 
neutrino 2/3 decay of ^^°Gd is strongly forbidden due to heavy deformation of this nucleus. 
Meanwhile, the suppression of the Oi/2/5 decay mode would be not so strong due to different 
sets of intermediate states involved in both transitions. Therefore the energy region of Oz/2/5 
signal of ^^°Gd could be free of the background from 2i/2/5 decays, which is very serious problem 
for 2f3 detectors with the poor energy resolution |0 . Thirdly, the natural abundance of ^^°Gd 
is rather large (21.86% |^) allowing to build up the high sensitive apparatus with natural 
GSO crystals. 

The present paper describes the new and further improved half-life limits on 2/5 decay of 
^^°Gd obtained with the help of 95 cm^ GSO crystal scintillator and with about 3 times larger 
running time than in ref. [^, where preliminary results of this experiment have been already 
published. 



2 Set up, background measurements and data analysis 

Cerium-doped gadolinium silicate crystal (5.4 cm long, 4.7 cm in diameter) grown by Czochral- 
ski method was used in the measurements. The mass of the crystal is 635 g, and the number 
of ^^°Gd nuclei is 3.951x10^^. The first 630 h of measurements had been carried out with the 
mass of the crystal equal to 698 g, then its side surface had been ground on 1-1.5 mm. 

The experiment was performed in the Solotvina Underground Laboratory (SUL) of the INR 
in a salt mine 430 m underground (^1000 m w. e., cosmic muon flux 1.7x10"^ cm~^ s~^, 
neutron flux <2.7xl0^® cm^^ s~^, radon concentration in air <30 Bq m~'^) ^3]. In the low 
background installation the GSO crystal is viewed by the photomultiplier FEU-110 through a 
plastic light-guide 8.6 cm in diameter and 18.2 cm long. The energy resolution of the detector 
was measured in the energy region 60 — 2615 keV by using 7 lines of ^^Na, ^^^Cs, ^°^Bi, ^^^Ra, 
^^^Th and ^^^Am sources. As an example, the resolution equals 16.8%, 13.5%, 11.2% and 10.7% 
at the energy 662, 1064, 1770 and 2615 keV, respectively. In course of measurement the energy 
calibration was carried out with ^''^Bi source weekly. The passive shield made of high purity 
copper (5 cm thickness), mercury (7 cm), and lead (15 cm) surrounds the GSO scintillator to 
reduce the external background. Event-by-event data acquisition system consists of the IBM 
PC compatible personal computer and CAMAC crate with electronic units, which allow to 



record the amplitude (energy) and arrival time of each event [p^ . 

^For example, the first interesting results were obtained recently for 2/3 decay processes in ^°Ca and ''^Ca 
with the help of newly developed low radioactive CaF2 (Eu) crystal scintillators IQ . 
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Total statistics collected in the experiment is 13949 h (1.015 yr kg of exposure). The 
measured background spectrum of the GSO crystal is depicted in fig. 1, where the following 
pecuharities exist: the clear peak at the energy 420 keV, the comparatively wide peak at the 
energy around 1050 keV and two broad distributions dropped down at the energies 2.4 and 
5.5 MeV. Taking into account the relative light yield for a particles as compared with that for 
electrons (a//5 ratio) for the GSO scintillator^, the first peak is attributed to a particles of ^^^Gd 
(Ti/2 = 1.08x10^^ yr; = 2140 keV; abundance 6 = 0.20%) and ^^^Sm (T1/2 = 1.06 xlO^^ yr; 
Ea = 2233 keV; 6 = 15%; samarium can be present as impurity of the GSO crystal at the level 
of ~8 ppm [0). The peak near 1050 keV as well as the broad distribution up to the energy 
2.4 MeV is mainly due to the radioactive contamination of the crystal by the nuclides from the 
^^^Th, 235u g^]^^ 238-[j fajxiilies. The second distribution up to the energy 5.5 MeV is caused by 
decays of ^^^Th daughter isotopes: a) (3 decay of ^osjl (Q^j = 5.00 MeV); b) (3 decay of ^i^Bi 
{Q/3 = 2.25 MeV) followed by a decay of its daughter ^^^Fo (T1/2 = 0.3 /xs; E„ = 8.78 MeV or 
^2.7 MeV in p scale). 
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Figure 1: Background spectrum of the GSO crystal (95 cm^) collected during 13949 h. In the 
insert: part of the measured distribution in the energy interval 600 - 3000 keV (open circles) 
together with the residual after subtraction (see text) of the decays from the ^^^Ra, ^^'^Ac and 
^^^Th intrinsic contamination of the crystal (filled circles) and fitting curve (solid line). 



To recognize and reduce the background from the intrinsic radioactive impurities of the 
crystal, the off-line analysis of the arrival times of measured events was fulfilled (as described 
in details elsewhere |^ )• Using this method fast sequences of decays belonging to the 
natural radioactive chains were searched for, as for example the sequence of two a decays from 

The energy dependence of the a/ (3 ratio was determined by using the a peaks of ^"Po, 2i5po^ 2i6pQ g^-,^^ 
^^°Rn from the internal contamination of the crystal as following: a/P = 0.152 + 0.01765 • Ea, where Ea is in 
MeV. The mentioned peaks were selected from the background with the help of the time-amplitude analysis 
described below. 
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the 232Th family: ^^^Rn {Q^ = 6.41 MeV, T1/2 = 55.6 s) ^ ^lep^ (g^ _ g ^eV, Ty2 = 0.145 
s) ^^^Pb. Because the energy of ^^°Rn a particles corresponds to 1.7 MeV in (3 scale of the 
GSO detector, the events within the energy region 1.2 — 2.2 MeV were used as triggers. Then 
all events following the triggers in the time interval 10 — 1000 ms (it contains the part rjt = 
0.945 of the total number of the ^^^Po decays) were selected. As an example, the spectra of the 
^^°Rn and ^^^Po a decays obtained in this way - as well as the distribution of the time intervals 
between the first and second events - are presented in fig. 2. It is evident from this figure 
that selected spectra and time distribution are in the excellent agreement with those expected 
from a particles of ^^°Rn and ^^^Po. Taking into account the efficiency of the time-amplitude 
analysis, the number of accidental coincidences, and the interfered chain of ^^^Rn — >^^^Po (it is 
selected by the applied procedure with 2.04% efficiency), the ^^^Th activity in the GSO crystal 
was determined as 2.287(13) mBq/kg. Then on the next step of analysis the fast couples found 
(^^°Rn and ^^^Po) were used as triggers to search for preceding a decays of ^^^Ra {Qa = 5.79 
MeV, Ti/2 = 3.66 d). The time window was set as 1 — 30 s (it contains rjt = 0.30 of ^^°Rn 
decays). Resulting distribution which includes the accidental coincidences (calculated ratio of 
the effect to the accidental background is equal to 0.825) is also in a good agreement with the 
expected a peak of ^^''Ra. 
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Figure 2: The energy spectra of the first and second a particles in the decay chain ^^°Rn 
^^^Po ^^^Pb which were found with by means of the time-amplitude analysis of the data 
recorded over 8609 h. In the insert: the distribution of the time interval between the first 
and second events together with its fit (solid line) by the sum of exponent (dashed line) with 
Ti/2 = 0.145 s (table value Ty2 = 0.145(2) s ^) and background (dotted line). 



The same technique was applied to the sequences of a decays from the ^^^U family: ^^^Ra 
{Qa = 5.98 MeV, T1/2 = 11.44 d) ^ ^^^Rn (g„ = 6.95 MeV, T1/2 = 3.96 s) ^ ^isp^ (q^ _ 7 53 
MeV, Ti/2 = 1.78 ms) ^"Pb. For the fast couple (^^^Rn and ^^^Po), events within 1.3 - 2.4 
MeV were used as triggers and the time interval of 0.5 — 10 ms (containing rjt = 0.803 of ^^^Po 
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decays with energy between 1.5 and 2.6 MeV) was chosen^. The obtained a peaks correspond 
to activity of 0.948(9) mBq/kg for the ^^''Ac impurity in the crystal Then the procedure 
analogous to that described for the ^^"^Ra was apphed to find the preceding ^^^Ra a decays. 

For the analysis of the ^^^Ra chain (^^^U family) the following sequence of /3 and a decays 
was used: ^i^Bi (Q^ = 3.27 MeV, T1/2 = 19.9 m) ^ ^Mp^ ^ 7 gg ^^y^ y^^^ ^ 3 ^ 

^^"Pb. For the first event the lower energy threshold was set at 0.5 MeV, while for the second 
decay the energy window 1.3 — 3.0 MeV was chosen. Time interval of 2 — 500 /xs ijit = 0.872 
of ^^^Po decays) was used. The selection efficiency is also decreased a little by the energy 
threshold applied to the first event (in the procedures described above the selection efficiency 
of the energy windows rjE is equal to 1). By the Monte Carlo simulation the part rjE of the 
^^'^Bi spectrum above 500 keV (as compared with the total spectrum) has been determined as 
rjE = 0.794. The obtained spectra for ^^'^BiQ and ^^^Po are shown in fig. 3 and lead to the ^^^Ra 
activity in the GSO crystal equal to 0.271(4) mBq/kg. 
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Figure 3: The energy spectra of the sequence of (3 and a decays in the decay chain ^^^Bi 
2i4po 2^opb which were found by means of the time-amplitude analysis of 8609 h data. In 
the insert: the distribution of the time interval between the first and second events together 
with its fit (solid line) by the sum of exponent (dashed line) with T1/2 = 129 /zs (table value 
Ti/2 = 164.3(20) /is [^) and exponent with T1/2 = 1.78 ms related with the chain ^^^Rn — > 
2i5po ^ 2iip^ (dotted line). 



Besides determination of the background components, the time-amplitude analysis is used 

^The used procedure selects also pairs ^^"Rn — > ^^^Po and ^^"^Bi — > ^^"'Po with efficiencies 4.64% and 12.08%, 
respectively. These contributions are taken into account in the calculation of activities. 

^Peak observed in the f3 spectrum of ^^^Bi (see fig. 3) is the part (17.7%) of ^^^Rn a decays (from ^35^ 
family), which corresponds to the chosen time interval 2—500 /iS. Due to known activity of ^^''Ac the contribution 
of ^^^Rn ^^^Po (as well as 0.24% of the chain ^^"Rn ^ ^^^Po) was calculated accurately and subtracted from 
the activities of ^^^Bi and ^^"'Po. 
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for reduction of the background. In this case, the attention is paid, first, to selection of time- 
correlated decays with minimal loss of the exposure, and, secondly, to prevention of the double 
selection of the events. With these aims four "removing" procedures were applied (the sub- 
scripts numerate events in the chain): 

(i) El > 500 keV; At = 2 - 500 /is; E2 = 1000 - 3000 keV. 

This cut removes mainly decays of ^^^Bi and ^^^Po (774 = 0.872 , rjE^ = 0.794, riE2 = 1.00, 
thus the total efficiency of selection is 77 = rjtrjE^'nE2 = 0.692), while efficiencies for the pairs of 
2^^Rn - 2i5po and ^^^Rn -^^^Po are rather low: 0.177 and 0.024, respectively. 

(ii) El, E2 = 1200 - 3000 keV; At = 500 /is - 1 s. 

The corresponding efficiencies are: ?7(^^^Rn - ^^^Po) = 0.823; //(^^^Rn - ^^^Po) = 0.989; and 
?7(^^^Bi - ^^^Po) = 0.067. It is also necessary to take into account the pairs of ^^^Ra - ^^^Rn 
{r] = 0.161) and ^s^Ra - 220j^^ ^ 0.0124). 

(iii) In the third cut a fast chain of decays (^^°Rn - ^^^Po) is used as a trigger to select the 
previous event of ^^^Ra. The parameters taken for ^^^Ra selection are the following: Ei >1000 
keV; Ati2 = 1 - 30 s; = 1200 - 2200 keV; Atss = 0.01 - 0.5 s; E^ = 1300 - 2400 keV. The 
final efficiency for selection of ^^^Ra is rj = 0.258. 

(iv) In the fourth cut procedure a fast chain of decays (^^^Rn - ^^^Po) is used as a trigger 
to select the previous event of ^^^Ra. The parameters are: Ei >1000 keV; Ati2 = 1 — 10 s; 
E2 = 1300 - 2400 keV; Atas = 500 /is -10 ms; E3 = 1500 - 2600 keV. It removes 53.4% (r/ = 
0.534) of 223Ra decays. 

Since the time windows of the cuts are not overlapping, the probability that one event could 
be selected more than once is rather low. Indeed, from the total number of 345169 events being 
attributed to the time-correlated background, only 318 (i. e. 0.092%) are double-selected. 

The part of the measured distribution in the energy interval 600 — 3000 keV is depicted in 
the insert of fig. 1 (open circles) together with the resulting spectrum after subtraction of the 
selected decays (filled circles). It corresponds to 0.969 yr-kg of exposure (~ 13400 h), or 95.7% 
of the initial value that was decreased by the "removing" procedures. At the same time, due to 
these procedures the background rate in the vicinity of the energy release of the Oz/2/3 decay of 
^^"Gd (1648 — 1856 keV) has been reduced by 2.3 times up to the value of 1.01(1) cpd/keV-kg. 



3 Background simulation 

With the aim to evaluate the 2/3 decay processes in ^^°Gd the measured background spectrum 
(after subtraction of the time-correlated decays) was simulated with the help of GEANT3.21 
package ||26|. The event generator DECAY4 was used to describe initial kinematics of 
decays (number and types of emitted particles, their energies, directions of movement and 
times of emission). It takes into account decays to ground state as well as to excited levels 



of daughter nuclei with the subsequent complex de-excitation process p5[|. The possibilities of 
emission of conversion electrons and e"'"e~ pairs instead of 7 quanta in nuclear transitions and 
the angular correlation between emitted particles are also taken into consideration. 

The background model was built up as a result of the procedure, in which the experimental 
spectrum was fitted by the sum of simulated response functions. The coefficients of the latest 
were determined on the basis of the following data: 

i) Activity values for ^^^Th, ^^'^Ac and ^^^Ra (and their short-lived daughters) that are 
present as intrinsic contamination in the GSO crystal. These activities were determined firmly 
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and accurately (with uncertainty less than 1%) with the help of the time-amplitude analysis, 
as described above. The part of decays of ^^^Ra, ^^^Rn, ^^^Fo; ^^^Ra, ^^^Rn, ^^^Fo; ^^^Bi 
and ^^^Po was removed from the background spectrum, and the remaining part is calculated 
with high precision for all these nuclides. Thus, one can describe the spectrum of, for example, 
^^®Th+ daughters as a sum of the simulated spectra with exactly known areas. If we suppose the 
secular equilibrium within the natural radioactive chains for these contaminations, the activities 
of the remaining long-lived members of the ^'^^Th, ^ss-y ^^^^ 238-u faj^iUes would be known too. 
However, it is known that chemical procedures and crystal growth usually break the equilibrium 
in the natural radioactive series. To account this possibility the activities of the mentioned 
remaining members of the radioactive chains (namely, ^^^Th, ^^^Ra, 235u_j_238-[j_j_234-[j^ ^^^Pa, 
^^°Th, and ^^°Pb) were taken as free parameters for the fitting procedure^. 

ii) The radioactive impurities of the photomultiplier (PMT), which are the main source of 
the external background. Their values for PMT FEU-110 were measured previously as 



3.0(3) Bq (^oK), 0.8(2) Bq (^^^Ra) and 0.17(7) Bq C^^Th). In the fit these activities were taken 
as free parameters varied within their errors. 

Besides, simulated spectra of "^^K and ^^^La - natural radionuclides which could be present 
in the GSO crystal - were included into the fitting procedure too. The last component of the 
background model is the exponential function (with two free parameters) which describes the 
residual external background (multiple scattering of 7 quanta, influence of weak neutron flux 
and so on). The exponential behavior of this component was confirmed by the measurements 
with the high radiopurity CdW04 crystal scintillator (454 g) performed in the same set up 



The fit of the experimental spectrum in the energy region 0.1 — 3.0 MeV by the sum of the 
described components gives the following activities of the additional intrinsic contamination of 
the GSO crystal: < 14 mBq/kg; ^^SLa < 55 mBq/kgQ ^^^Th < 6.5 mBq/kg; Ra < 9 
mBq/kg; ^ssy < 2 mBq/kg; ^^^Fa < 0.08 mBq/kg; ^^^Th < 9 mBq/kg; ^lOpb < 0.8 mBq/kg. 
The parameters of the exponent were also found in the fit procedure and its contribution to the 
experimental distribution was revealed to be small (~ 2% for the energy interval 1 — 2 MeV). 
For illustration the fitting curve (in the energy region 0.7 — 2.4 MeV) is presented in the insert 
of fig. 1 together with the experimental data. 



4 Half-life limits on the 2f3 decay of ^^"Gd 

Since in the measured spectrum the Oz/2/3 decay peak of ^^°Gd is evidently absent, only the limit 
for the probability of this process can be set on the base of the experimental data. To estimate 
the half-life limit T1/2, we use the formula lim T1/2 = ln2 ■ rj ■ N ■ t/\imS, where r] is the 
detection efficiency, N is the number of ^^°Gd nuclei, t is the measuring time and lim S is the 
number of effect's events which can be excluded with a given confidence level. To calculate the 
values of t] and lim 5, the response function of the GSO detector for the effect being sought has 
been simulated with the help of GEANT3.21 and DECAY4 programs. It was found that for the 
Oz/2/5 decay the response function is the Gaussian centered at 1730 keV and its width (FWHM) 
equals 176 keV. The edge effects (escape of one or both electrons and bremsstrahlung quanta 

^Three isotopes of uranium ^^^U, ^ssy g^^j ^34^ j^q^ chemically separable, and their relative activities 
(0.046:1:1) do not change, so we can take only one parameter for these nuclides. 

^This limit of ^"^^La activity corresponds to the possible La impurity in the GSO crystal at the level of 67 
ppm, which does not contradict with the results of the chemical analysis. 
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from the crystal) remove from the peak ^ 5% of events, thus t] = 0.95. The hmS* values were 
determined in two ways. First, by using the so called "one a approach", in which the excluded 
number of effect's events is estimated simply as square root of the number of background counts 
in a suitably chosen energy window AE. Notwithstanding its simplicity this method gives the 
right scale of the sensitivity of the experiment. For instance, in the measured spectrum within 
the energy interval 1648 — 1856 keV (it contains 82% of the expected peak area) there are 
74500 counts; thus, the square root estimate gives limS" = 273 events. Using this value of 
limS" , the total exposure related to ^^°Gd nuclei {N ■ t = 6.04 x 10^^ nuclei yr), and the 
calculated efficiency [r] = 0.78), we obtain the half-life limit: T1/2 > 1.2 x 10^^ yr (68% C.L.). 
Further, the lim S value was determined by using the standard least squares procedure, where 
the experimental energy distribution in the vicinity of the peak searched for was fitted by the 
sum of the background model (as described above) and effect's peak being sought. It should 
be stressed that for the energy interval of interest the main important contributions (~73% of 
the experimental spectrum within 1600 - 1900 keV) are the activities of ^^^Ra, ^^^Ac and ^^^Th 
(and their short-lived daughters) from the intrinsic contamination of the GSO crystal, whose 
values were determined accurately (better than 1% uncertainty). 

As a result of the fit procedure in the energy region 1.3-2.1 MeV, the obtained area for 
the Oz/2/3 decay peak is —160 ± 233 counts (x^ value equals to 0.85), thus giving no evidence 
for the effect. The difference between measured and simulated spectra in the energy region of 
the hypothetical 0t^2/? decay peak of ^^''Gd is shown in fig. 4, where the excluded effect is also 
depicted (solid line). The number of effect's events, which can be excluded with 90%(68%) 
C.L. was calculated |2| as 298(169). It gives the half-life limit: 
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Figure 4: The residual between the experimental and simulated spectra in the energy region of 
the 01^2(3 decay of ^^°Gd (open circles with error bars). The solid line represents the excluded 
Oz/2/5 peak with T^/^ = 1-3 x lO^i yr (90% C.L.). 

Ti/2(0z/2/5) > 1.3(2.3) X 10^1 yr at 90%(68%) C.L. 
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Comparing our limit with the theoretical calculations [IS], one can compute the restriction 
on the neutrino mass (m^) < 26(19) eV at 90%(68%) C.L. 

The same T1/2 limit has been set for the transition to the first excited level of ^^°Dy (2+, 
87 keV) because 87 keV 7 quanta following this process will be almost fully absorbed inside 
the scintillator. The obtained limits are several times larger than those already published 



15, 16, 17, 22 



The two neutrino 2(3 decay rate was evaluated in two ways. For the very conservative 
estimate the model spectra of the exactly measured contamination in the crystal and PMT were 
subtracted from the experimental data. The residual in the chosen energy region was equated 
to the simulated 2v2(3 decay distribution and the latest was taken as an excluded effect with 
lim S = 7.6x 10^ events corresponding to the Ty2{S-^-) > 5.5x 10^^ yr (99% C.L.) for 2i/2/5 decay 
of ^^°Gd. However, within this simple approach it was impossible to reproduce adequately the 
measured spectrum, thus the background model and the fitting procedure described above were 
applied for the 2^ decay mode too. The set of fits were performed by changing the fitting energy 
region from (100 — 760) keV to (2400 - 3000) keV. The maximum value of an excluded effect 
(lim 5* = 2.2x10^ events at 90% C.L.) was found for the energy interval 760 - 2600 keV. The 
corresponding Monte Carlo simulated spectrum of 2z/2/3 decay of ^^°Gd (g.s.— >g.s.) is shown 
in fig. 5 together with the fitting curve and the most important background components. It is 
visible from this figure that our background model reproduces the experimental data quite well 
even outside the region of the fit. The final lower limit for the process searched for is equal to: 
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Figure 5: The background spectrum of the GSO detector for 0.969 yr-kg of exposure (points) 
and the model of background (solid line) obtained by the fitting procedure in 760 — 2600 keV 
energy interval (see text). The most important internal (^°K and sum of ^^^U, ^^^U, ^^^Th) and 
external (7 radiation from PMT) components of background are shown. The excluded with 
90% C.L. distribution of 21/2/3 decay of ^^°Gd to the ground level of ^^°Dy corresponds to the 
half-life limit T^Ug-s.) = 1.9 x 10^^ yr. 
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Tf/2(g-s.) > 1.9(3.1) X 10^9 yr at 90%(68%) C.L. 

The same method gives the hmit for 2z/2/? transition to the first (2"'") excited level of ^^°Dy0: 

7?;2(2+) > 2.1(3.4) X 10^^ yr at 90%(68%) C.L. 

The similar fitting procedure was used to set limits on double beta decays with one or two 
Majoron emission. The estimated restrictions on half- life are: 

T°/7 > 3.5(5.3) X lO^s yr at 90%(68%) C.L. 
T°7^ > 1.3(2.0) X 10^9 yr at 90%(68%) C.L. 

5 Limits on 2(3 decay processes of Ce isotopes 

The concentration of cerium in the GSO(Ce) crystal (0.8%) is known from the crystal growth 
conditions and from the results of the chemical analysis. It allows to search for the 2(3 processes 
in three cerium isotopes: double positron decay (2/5+), or electron capture and positron decay 
(e/3+), or double electron capture (2e) in ^^^Ce (mass difference between parent and daughter 
atoms AMyi = 2397(48) keV; abundance of parent nuclide 5 = 0.185%); double electron capture 
in i38Ce {AMa = 693(11) keV; S = 0.251% ); and 2/?" decay in ^^^Ce (AMa = 1417(2) keV; 

6 = 11.114%). The total numbers of ^^^Ce, ^^^Ce and ^^^Ce nuclei in the crystal are 4.1x10^^, 
5.4x10^^ and 2.4x10^^, respectively. The response functions of the detector for the different 
possible 2/3 processes in these cerium isotopes were simulated with the help of GEANT3.21 
code and event generator DECAY4. As an example, the simulated energy distributions for 
neutrinoless and two-neutrino 2/3^ (K/5"'"and 2K) decays of ^^^Ce are presented in fig. 6. Limits 
on half-lives with respect to different modes of 2/3 processes of cerium isotopes were calculated 
in the way analogous to that for ^^'^Gd - they are summarized in the Table. Except two limits 
for 0z/2/3+ decay of ^^^Ce and Oz/2/5~ decay of ^^^Ce (already obtained by using CeFs scintillators 



30|), all presented results for cerium isotopes are set for the first time. 



6 Conclusions 

The present experiment performed in the Solotvina Underground Laboratory with the help of 
GSO scintillator has shown "pro et contra" of these crystals for the advanced double f3 decay 
study of ^''"Gd. 

The current level of sensitivity is limited mainly by the intrinsic contamination of the GSO 
crystal (^^^Th, ^ss-jj 238-[j families), therefore further development (which is in progress 
now) to remove these impurities from the crystals is at most importance. Careful purification 
of raw materials from actinides and their daughters (technically available now) could decrease 
radioactive contaminations of the GSO crystals by two-three orders of magnitude down to 
the level of several /xBq/kg^. Another possibility to reduce background is the pulse-shape 

®It should be noted that shape of distribution for 2^2(3 decay to 2+ level differs from 2i^2/3 transition to 
ground (0+) state [|. 

^Such a radiopurity has been already reached in the CdW04 crystal scintillators used for the 2(3 decay study 

of ii^cd niini. 
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Figure 6: The response functions of the GSO(Ce) detector for neutrinoless and two-neutrino 
2(3+ (K/5+ and 2K) decays of ^^^Ce simulated with the help of GEANT3.21 code and event 
generator DECAY4. 



discrimination. Our preliminary measurement shows the difference of decay times for pulses 
induced in GSO by a particles and 7 quanta: 50.8(4) and 48.3(4) ns, respectively]^ Although 
this difference is not so large, the proper pulse shape discrimination technique, which is under 
development now, would allow to eliminate a background from decays of ^^^Sm, ^^^Gd and 
members of ^^^Th, ^^^U and ^^'^U families. 

On the other hand, due to the high abundance of ^^°Gd (~ 22%) the GSO crystals could be 
grown up from natural Gd, thus such detectors should be well less expensive than those made 
of enriched 2(3 decay candidate isotopes. Therefore, the enforcement of the large scale and 
high sensitivity experiment with ^^°Gd would be really possible by using the GSO multi-crystal 
array with total mass of about 2000 kg (;^400 kg of ^^°Gd). 

The background of such a detector could be reduced further by placing these crystals into a 
high purity liquid (water or scintillator)!^ serving as shield and light guide simultaneously. The 
similar idea has been recently proposed in CAMEO project with aim to study 2(3 decay of 



^^^Cd by using ^^^CdW04 crystals placed in the liquid scintillator of the BOREXINO Counting 
Test Facility (CTF). It was evidently demonstrated by pilot measurements with ^^^Cd and 
by Monte Carlo simulation that sensitivity of the CAMEO experiment (in term of the T1/2 
limit for 0^2(3 decay) with f» 100 kg of ^^^CdW04 crystals is ^10^^ yr which translates to the 
constraint on the neutrino mass m^, < 0.06 eV [0. Moreover, the strong dependence of the 

These results are in a good agreement with ref. ||3^ . 
-'^^The existing and future large underground neutrino detectors (SNO |Q , BOREXINO KamLand js^) 
could be appropriate for our proposal with one-two tons of the GSO crystals. The latest located in the water 
or liquid scintillator would be homogeneously spread out on a sphere with diameter 3-4 m and viewed by the 
distant PMTs. 
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light collected by each PMT versus coordinate of the emitting source in the crystal has been 
found. Such a dependence is explained by difference of the refraction indexes of crystal [n = 
2.3 for CdW04) and liquid scintillator {n = 1.58), which leads to the redistribution between 
reflected and refracted light due to change of the source position. By means of the GEANT 
Monte Carlo simulation it was shown that spatial resolution of 1- 5 mm (depending on the event 
location and the energy deposit) can be reached with CdW04 crystals (07 x 9 cm) placed in the 
liquid scintillator of the CTF and viewed by 200 distant PMTs With the GSO detectors 



(refractive index n = 1.85) placed in a liquid (n ^ 1.5), the simulation of light propagation 
gives the spatial resolution in the range of 4 - 10 mm. Anyhow, it would certainly allow to 
reduce background in the energy region of interest additionally (roughly by factor of 10 - 50). 

We estimate that sensitivity of the experiment with about two tons of the GSO crystals 
(placed in the SNO or BOREXINO set ups) and for 5-10 years of exposition would be of 
the order of ~2xl0^^ yr, hence the restriction on the Majorana neutrino mass can be 
reduced down to rriiy < 0.07 eV. It is comparable with the sensitivities of the recently proposed 
large scale projects for 2/3 decay study, like MOON |0[, EXO gll, CAMEO fS^, CUORE @; 



MAJORANA g^, GENIUS whose results could provide crucial tests of the certain key 
problems and theoretical models of the modern astroparticle physics. Note, however, that cost 
of the GSO experiment would be well lower than those of mentioned projects. 

The authors would like to thank A.Sh. Georgadze, B.N. Kropivyansky, A.S. Nikolaiko, and 
S.Yu. Zdesenko for their participation in the measurements. 
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Table 1: Half-life limits on the 2(3 decay processes in ^^OGd, ^^^Ce, i^«Ce and ^^^Ce 
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138 ( 



Nuclide 
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present work other works, (CL. 
90(68)% CL. 
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